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Abstract 

Accurate "first-principle" expressions for the excess free energy F^y^ and in- 
ternal energy Ue^ of the classical one-component plasma (OCP) are obtained. 
We use the Hubbard-Schofield transformation that maps the OCP Hamilto- 
nian onto the Ising-like Hamiltonian, with coefficients expressed in terms of 
equilibrium correlation functions of a reference system. We use the ideal gas 
as a reference system for which all the correlation functions are known. Ex- 
plicit calculations are performed with the high-order terms in the Ising-like 
Hamiltonian omitted. For small values of the plasma parameter T the Debye- 
Huckel result for F^x and f/gx is recovered. For F » 1, these depend linearly 
on F in accordance with the Monte Carlo findings for the OCP. The MC data 
for the internal energy are reproduced fairly well by the obtained analytical 
expression. 

I. INTRODUCTION 



The one component plasma (OCP) model is one of the basic models in the condensed 
matter physics Besides of its direct astrophysical applications p|, to model the ionized 

1 



matter in white dwarfs , outer layers of neutron stars and interiors of heavy planets , 
OCP is widely used as a reference model for a variety of systems, ranging from alkali metals 
to colloidal solutions P-pJ]]. The application of the OCP is not restricted by classical 
systems; it is also used when quantum effects are important |]12| . 

The OCP model is formulated as a system of point particles, interacting via the Coulom- 
bic potential, which move in a uniform neutralizing background ll|J^. All thermodynamic 
properties of the OCP depend only on the dimensionless plasma parameter F = /b/^c, where 
Ib = e^/kBT is the Bjerrum length (e is the charge of the particles, is the Boltzmann's 
constant, T is the temperature) and Oc = (3/47rp)^/'^ is the ion-sphere radius with p = N/il 
being the concentration of particles (A^ is the number of particles, Q is the volume of the 
system) . 

For small values of F, which correspond to a hot and/or dilute system, the Debye-Huckel 
theory accurately describes the thermodynamic properties of the OCP. For the excess free 
energy and internal energy this gives in the F ^ limit: 



ksTN V3 ' keTN 2 

Abe expansion |T3] for the classical OCP provides next terms for that limit (in the present 
study we consider the classical case): 

^'"^ -^T''^ - "-T' - (31ogF - 1) + • ■ • (2) 
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where c = f log3 + §7 - 1 = 1.101762 . . ., and 7 is the Euler's constant. The analytical Abe 



expansion seems to be fairly accurate for F up to 0.1 |14|. The next few terms for the F — > 
expansion calculated in Ref. |T3| allow to use the small-F expansion up to F < 0.4. 

For larger values of F the OCP was studied numerically, by means of integral equa- 
tions, such as Percus-Yevick, hypernetted-chain equations [|rB|, and (most successively) by 
modified hypernetted-chain equation . Extensive numerical studies have been performed 
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by Monte Carlo (MC) P,p!^,p!8|-pT| and Molecular Dynamics p2| technique. To fit avail- 



able "experimental" data for the excess thermodynamic functions simple analytical fits were 
proposed ||23|| : 

^'^'^ AT + Br + C (4) 
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AT + -r + {3 + C)\ogT-D (5) 
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with A = -0.8992, B = 0.596, s = 0.3253, C = -0.268 and D = A + B/s + 1.1516. Eqs.® 
and d^) are fairly accurate for the interval 1 < F < 220 but unfortunately they do 



not give a correctly behavior at F — ^> 0. Fade approximants for the t/ex suggested in p4|j2l 
remedy the failure of (HI) at small F; one thus obtains very precise description for the whole 



interval < F < 200 p5 



To describe the limit of very large F a perturbation theory was proposed p6[; it agrees well 
with the MC data for large F, but also lack the proper small-F behavior. The correct Debye- 
Huckel behavior at F ^ together with a reasonable ~ 10% accuracy for < F < 100 has 
been obtained in a simple semiphenomenological "Debye-Huckel plus hole" (DHH) theory 
| |27| |; it fails, however, for F > 125 A modified DHH theory proposed recently ||2S| using 
only one fitting parameter accurately reproduces the MC results for the large values of F 
(1 < F < 200) and demonstrates a correct behavior at F 0. 

Thus, up to now, no "first-principle" theory of OCP exists which describes accurately the 
thermodynamic properties in the whole range of F from the Debye-Huckel limit F — > up 
to F ^ 1 limit, where the Wigner crystallization occurs. In the present study we report 
a simple "first-principle" equation of state for the OCP which has the correct Debye-Huckel 
behavior for small F and demonstrates a linear dependence on F for F ^ 1. It reproduces 
within 1 — 2% accuracy the experimental data for the most of range of F (0 < F < 250) and 
has a typical deviation of the order of 2 — 5% (with a maximal one ~ 8%) for 0.01 < F < 1. 
The rest of the paper is organized as follows: in the next Sec. II we consider the Hubbard- 
Schofield transformation that maps the OCP Hamiltonian onto the Ising-like Hamiltonian 
and calculate the coefficients of the effective Hamiltonian. In this section we also present the 



field-theoretical formulation for the statistical sum of the OCP which directly follows from 
the transformation used. Within the Gaussian approximation for the effective Hamiltonian 
we derive the equation of state for the OCP. In Sec. Ill we discuss the equation of state 
obtained and compare the analytical results for the internal excess free energy with the 
available Monte Carlo data for the OCP. In the last Sec. IV we summarize our findings. 



II. EFFECTIVE HAMILTONIAN AND EQUATION OF STATE FOR THE OCP 

We start from the OCP Hamiltonian which may be written as follows {l3~^ = ksT): 

H = l/3-'Y.'^k (PkP-k - p) + ^^id (6) 
^ k 

where the first term in the right-hand side of Eq.(^) refers to the Coulombic interactions, 
written in terms of the collective variables, 

Pk = ^Ee-^''^^ (7) 

V " j = l 

where denotes coordinate of j-th particle, z/^ = AttIb / k"^ and Hi^ is the ideal-gas part of the 
Hamiltonian. Summation in Eq.(^ is to be performed over the wave- vectors k = {kx-, ky, kz} 
with ki = {i = X, y, z), where k are integers, = Q, and the prime over the sum denotes 



that the term with k = is excluded 29 . 



A. Hubbard- Schofield transformation 



The configurational integral may be then written in terms of the configurational integral 



of the reference (ideal gas) system Q ^ pO| , |3T 



as 



Q = ^exp|-i^Vfc(pkP-k-p)|^ Qr 



(8) 



where ((■ ■ = / dr'^ (■ ■ ■) denotes the averaging over the reference system. In accor- 
dance with the Hubbard-Schofield scheme ||3^ we use the identity. 
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1 r+oo 

exp(--a^a;^) = (27ra^)"^/^ / exp(--?/V«^ + ixy)dy 

Z J — oo Z 

and arrive after some algebra at: 

Q = QrJ n'ckC?</'kexp|-i^'z/^VkV'-k| ^exp |i^'pk</'-k| 



(9) 



R 

1/2 



where Ck = {2TiVk)~ ea'^'-'^, and where the integration is to be performed under the restric- 
tion, = ipl. [ifl^ is the complex conjugate of <y9k) HI]- Applying the cumulant theorem 
to the factor (exp I]k'PkV^-k})ij one obtains: 

Q = Qr n'^kC^V^kC"^, with 
k 

oo ^ 

H=J2^^''^ J2 M„(ki,...k„)v5ki---V5k„ 

n=2 ki,...k„ 

U2 (ki, ka) = - (5ki+k2,o I + (pkiP-ki)cfl I 

u„ (ki, . . .k„) = -i" — p (pki • • •Pk„)cR ^>2 (10) 



here denotes cumulant average ||33| for the reference ideal-gas system. As it follows 



from Eqs.(|l^), Q is written in a similar way as the partition function for the magnetic 



system having the Ising-like Hamiltonian |^4|, where <y9^ are the Fourier components of the 



"spin-field", ip{r). Note that the coefficients of the effective Hamiltonian (|T^) are expressed 
in terms of the correlation functions of the ideal-gas system, which are perfectly known (e.g. 

I): 



^ (Pki • • -PkJeiJ = P'^ki + ...k„,0 (11) 

This yields the effective Hamiltonian: 

1 / r -1 °° i"fi^~t / 

^ = 2 ^ r + (47r/B)-^/c2J ipi,ip_k - j — Pllki^^k^^i^i ■ ■ ■ <^k„5ki+...k„,o 



k n=3 

= nG + n, (12) 

where we write explicitly the Gaussian part Ti^: 



(13) 



Since all the coefficients of the effective Hamiltonian are known, one can develop the usual 
perturbation expansion, with Tic being a reference part of the Hamiltonian and with Tii 



being perturbation (e.g. P^j36| j3 



A simple structure of the effective Hamiltonian 7i, Eq.(12) suggests a simple, closed- 
form field-theoretical formulation for the partition sum of the OCP, akin the Sine-Gordon 
representation of the Coulombic gas [3^, or the restricted primitive model of electrolytes 
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B. Field-theoretical model for the OCP 



Using the space-dependent field ip{v) 

ifijc) -- 



It is easy to show that under this transformation the terms in Eq.(p!2|) containing prod- 
ucts <^ki ■ ■ 'V'k^ give rise to the terms / (iry9"(r), the terms k'^'Pkf-k give rise to the term 
/ dr (Vif)'^, so that one can write 

'1 



n = p j dr 



oo -n 



;(47rM-MVv.)^-EV 

n=2 



(14) 



and recognize the expansion of e^'^ . 

Integration over the Fourier-components y^k in Eq.(|l^) converts into "field" -integration 
over the field '■p{r). It may be shown that the Jacobian of this transformation does not 
depend on (y9(r) and appears as a normalization constant. 

Noticing that 
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n'(27rz/fc) ^= /n^kexpl-ij^j^^^y^kV' 



(with the restriction ifl. = f-k) and that 



2kBT 



where U{0) is the so-called "self-energy" ||38| . |40[] , one arrives after some algebra at the field- 
theoretical expression of the statistical sum Zqcp of the OCP: 



-PF _ 7 



OCP 



/ exp {- / dv [i«:^2(Vy.)2 - } 
where Vtp denotes the "field" -integration, 



(15) 



n= dr 



(Vv?)^ + i(p — e 



tip 



(16) 



= 4:7iIbP = 47re^/9p is the inverse Debye screening length and 



/i = /iid + ksT - f/(0)/2. 



(17) 



Here PfXid = log(A'^/fi) is the ideal-gas chemical potential (A = h/ is the 
thermal wavelength). Deriving Eq.(|T5|) we write the statistical sum of the ideal gas as 

= exp [— log(A'^/r2)] = exp {— /3/iidp / rfr} and rescale the length pdr —>■ dr. 

As it follows from Eq. (|I5]) the potential function of the effective field-theoretical Hamilto- 
nian for the OCP reads V{!f) = i(p — e**^. This may be compared with the potential function 
Ks'g(v^) = cosy) of the Sine-Gordon model for the Coulombic gas [^. Note, that all imag- 
inary terms in Eq . ([TsD vanish after the field integration due to the symmetry properties of 
the Hamiltonian under the transformation — > —(p. 

Consider now the "saddle-point" approximation to the numerator in Eq. ([T5|) . The equa- 
tion for the "extremal" field which minimizes the effective Hamiltonian reads: 



4( 



I — le 



vp 



(18) 



Under transformation ip = ie^cp/kBT Eq.(|T8|) converts into 



47rp fe-^'^/'^s^ - 1 



(19) 



which is the usual mean-/ie/d Poisson-Boltzmann equation for the OCP. This is not surprising 
since the "saddle-point" approximation is essentially the mean-field one. 
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C. Equation of state for the OCP 



Now we concentrate on the Gaussian part of the effective Hamiltonian and show that 
even neglecting the non-Gaussian contribution to the effective Hamiltonian, one can obtain 
fairly accurate equation of state for the OCP, provided that a correct value of the ultraviolet 
cutoff in the k-space is used. The Gaussian approximation to Ti. corresponds actually to 
the Random Phase, or Debye-Hiickel approximation, (e.g. |^). Using He, Eq.(p!^ and 
performing (Gaussian) integration (e.g. over (/^k in (|10|), one easily finds for the excess 
free energy of the OCP: 



- /?Fex = logiQ/Qu) = 7;Y1' [P^k - log (1 + pz/^)] 

^ k 



(20) 



We argue that the summation in Eq.(^) should be carried out over a finite number of the 
wave- vectors k. In this we follow the Debye theory of the specific heat of solids (e.g. [0). 
Namely, we assume that the total number of degrees of freedom in the system, 3A^, should 
be equal to the total number of physically different modes with the wave- vectors k within 
the spherical shell of radius ko in the k-space. The number of modes is twice the number 
of the wave-vectors, since for each k one has a sine and cosine mode (the amplitude of the 
k-th mode is a complex number) |Q. Thus we obtain: 



2 — -An 



k^dk = 3N 



(21) 



where the factor Q/Stt^ appears when the integration in k-space is used instead of summa- 

1 /3 

tion. From Eq . (pTj) follows that k^ = (9p7r^) . A similar Debye-like scheme to find the 
cutoff ko was first proposed for plasma in , where somewhat different value of the cutoff 
wave-vector was reported. Using the ko obtained we write: 



N 



I Q 4:71 fk 



2 87r3 Jo 



k'dk 



log 1 + 



A;2 



K 



D 



k^ 









/ x'^dx 






Jo 




x"^ 



(22) 



where b 



1/3 



The last integral is easily calculated to obtain for the free energy 
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^ [log (1 + hV) - 6r] - - (6r)5 arctan f ^ ) (23) 



kBTN 4 ^ ° ^ ' ' 2' ' V V6r 

and for the internal energy: 



= J- (6r)i arctan f ^ ) (24) 

of the OCP. To obtain Eg. (^41) we use the relation f/cx = T dFc^/dT. Again we note that the 
same functional dependence for the excess internal energy (but with different coefficients) 
has been obtained in [44]. 

III. RESULTS AND DISCUSSION 

As it follows from Eqs.(|23|) and (24), for F — the Debye-Huckel behavior is recovered. 



On the other hand in the opposite limit F ^ 1 Eqs. (P3|) and (24) demonstrate the linear 
behavior on F with the leading term —AT in accordance with fits (^ and (^ to the MC 
data. The constant A reads 



q / 9 \ 



which is fairly close to the constant A = 0.899 ... of the fits (^), (^). 

In Fig.l the excess internal energy given by Eq.(^) is compared to the Monte Carlo data, 
taken from Ref. MM for 0.1 < F < 1 and Ref. EH^ for F > 1. Fig.2 shows the relative 



error of the analytical expression (0). As it follows from Fig.l and Fig.2 the equation of 
state is fairly accurate in the most of range of the plasma parameter. The maximal deviation 
of the analytical expression from the numerical data occurs at the intermediate values of the 
plasma parameter, 0.01 < F < 1. 

To analyze the reason of the enhanced deviation of the theoretical results from the 
numerical data at 0.1 < F < 0.5, one can address the small-F expansion of Uex [P3]- It 



was observed |^8|| that in spite of the correct Debye-Huckel limit, this does not reproduce 
correctly the next-order terms of the Abe expansion (^. This occurs due to limitations 
of the Gaussian approximation for the effective Hamiltonian. Since all the coefficients of 
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the effective Hamiltonian are Icnown, one can go beyond the Gaussian approximation and 
develop a usual perturbation scheme, based on the Gaussian Hamiltonian. 




FIG. 1. Shows the dependence of the excess internal energy of the OOP Uc^/NkBT on the 
plasma parameter F = Is/ac {Ib = /ksT, = 47r/9/3). Points give the Monte Carlo data 
(Ref. [14,18] for 0.1 < F < 1 and Ref. [6,20] for F > 1); lines represent the theoretical result, 
Eq.( 
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FIG. 2. Shows the relative error (%) of the analytical expression ( p^ for the excess internal 
energy of the OCP Uc^/NkBT as a function of the plasma parameter F. The Monte Carlo data 
of Ref. [14,18] for 0.1 < F < 1 and of Ref. [6,20] for F > 1 are used. For F < 0.1 Eq.(|2|) is 
compared with the Abe expansion Eq.(^) 
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In particular the equation of state in a form of the virial expansion may be recovered pT 



Unfortunately this expansion does not provide the closed analytical expression for the excess 
thermodynamic functions, which may be used with acceptable accuracy for all the range of 

r. 

Thus, dealing with a problem where the relevant range of plasma parameter is not 
known in advance one should preferably use the simple closed-form equations suggested in 
the present study. 



IV. CONCLUSION 

A "first-principles" equation of state for the one-component plasma is derived that has a 
correct Debye-Huckel behavior at the limit of small plasma parameter F and demonstrates a 
linear dependence on F at F ^ 1. The obtained coefficient 0.881 at the linear leading term 
is close to the corresponding coefficient 0.899 found in the Monte Carlo simulations. The 
simple analytical expression for the excess internal energy reproduces the MC data within 
1 — 2% accuracy for the most of range of F (0 < F < 250) and has a typical deviation of the 
order of 2 — 5% (with a maximal one ~ 8%) for 0.01 < F < 1. 
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